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Abstract 



In the large iV c (number of colors) limit, quenched QCD and QCD are iden- 
£S) ■ tical. This implies that, in the effective field theory framework, some of the 

low energy constants in (N c = 3) quenched QCD and QCD are the same up 
to higher-order corrections in the 1/N C expansion. Thus the calculation of 
the nonleptonic kaon decays relevant for the AI = 1/2 rule in the quenched 
approximation is expected to differ from the unquenched one by an 0(1/N C ) 
CN ■ correction. However, the calculation relevant to the CP-violation parame- 

ter e'/e would have a relatively big higher-order correction due to the large 
cancellation in the leading order. Some important weak matrix elements are 
poorly known that even constraints with 100% errors are interesting. In those 
cases, quenched calculations will be very useful. 
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Quantum chromodynamics (QCD) is the underlying theory of strong interaction. At 
high energies (3> 1 GeV), the consequences of QCD can be studied in systematic, per- 
turbative expansions. Good agreement with experiments is found in the electron-hadron 
deep inelastic scattering and hadron-hadron inelastic scattering (Drell-Yan process) [1]]. At 
low energies (< 1 GeV), however, QCD becomes non-perturbative. Quantitative results not 
based on symmetry properties along can only be calculated directly from QCD via computer 
simulations on a Euclidean lattice. 

In lattice QCD simulations, it is common to perform the quenched approximation which 
is to drop the internal quark loop contributions by setting the fermion determinant aris- 
ing from the path integral to be one. This approximation cuts down the computing time 
tremendously and by construction, the results obtained are close to the unquenched ones 
when the quark masses are heavy (^> Aqcd ~ 200 MeV) so that the internal quark loops 
are suppressed. Thus the quenched approximation is an efficient and reliable tool to use as 
long as all the quark masses are heavy. However, in the real world, the light quark (u, d 
and s) masses are less than Aq C d so the justification for quenching does not apply. In this 
range of quark mass, quenching effects can be systematically studied by quenched chiral 
perturbation theory (Q;\PT) @, an effective field theory of quenched QCD (QQCD). In 
QxPT, quenching could induce a double pole to the flavor singlet meson propagator and 
could make quenched quantities more singular than the unquenched ones in the infrared 
region. In the ultraviolet region, QCD and QQCD are different as well because QQCD can 
be considered as having infinite quark masses for the internal loops so its ultraviolet behavior 
is different from that of QCD. Those uncontrolled effects make it hard to see what we can 
learn about QCD from QQCD. 

However, there is a known connection between QCD and QQCD in the limit that the 
number of colors (N c ) becomes infinite; the large N c limit of QQCD is indistinguishable to the 
large N c limit of QCD. In this paper, we study the finite N c quenching errors systematically 
using effective field theories with 1/N C expansions. Generally, in the real world where N c = 3, 
some of the contributions which are formally higher order in the 1/N C expansion could 
become numerically important for certain observables near the chiral limit (zero quark mass 
limit). Fortunately, these big contributions are typically non-analytic in quark masses and 
can be computed in effective field theories of QCD and QQCD. After these big non-analytic 
contributions in QQCD are replaced by those in QCD, it is plausible that the errors which 
are formally 0(1/N C ) in the quenched approximation could imply numerically ~ 30% errors. 
Sometimes even smaller errors can be achieved by forming ratios of observables, such that 
the 0(1/N C ) errors from the leading terms in the chiral expansion cancel. The ratio of B 
meson decay constants /b s //b [@>F| is an explicit example. We then study the quenched 
calculations of the non-leptonic kaon decays K — ► tttt |5|||. We find that the quenching 
error in the part relevant to the AI = 1/2 rule is consistent with 0(1/N C ) ~ 30%. However, 
the CP violation parameter e'/e calculation could have a big quenching error due to the 
large cancellation in the leading terms in the 1/N C expansion. Finally, for illustration, we 
identify some cases that quenched calculations can be very useful. Those are calculations of 
important but poorly known weak matrix elements for which even 100% error constraints 
are interesting. 

It is known that the quenched world has some similarity to the large N c world. Take a 
pion two-point function (0 |7r(x)7r(0)| 0), for example: there are contributions from diagrams 
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FIG. 1. Feynman diagrams contributing to the pion two point function. The solid lines and 
curly lines denote the quark and gluon propagators respectively. 

shown in Fig. 1. For convenience, the quarks connected to the external sources are called 
"valence quarks" ; the quarks not connected to the external sources are called "sea quarks" . 
In the quenched approximation, the diagrams with sea quark loops such as the one in Fig. 
1(b) are zero. (Then the gluon-quark coupling in surviving diagrams such as the ones in 
Fig. l(a,c) are usually rescaled to mimic the situation that the sea quark masses are set to 
infinity. We will not do this rescaling at this moment but leave it for later discussion.) In 
contrast, in the large N c limit proposed by 't Hooft @, 



where g is the quark gluon coupling, both the Fig. 1(b) diagram with one sea quark loop 
and the "non-planar" diagram of Fig. 1(c) with one crossing in gluon lines are suppressed by 
one power of 1/N C compared to that of Fig. 1(a). The diagrams with more sea quark loops 
or more gluon crossings are further suppressed by even higher powers of l/N c . So, in the 
large N c limit of QCD, only the planar diagrams such as Fig.l(a) survive, while in QQCD 
both Figs. 1(a) and (c) survive. Now if we further take the large N c limit of QQCD, the non- 
planar diagrams vanish as well. From this we conclude that QQCD and QCD have the same 
large N c limit — at least in the mesonic two-point function case. It is easy to generalize the 
above observation to an arbitrary n-point function even involving baryons. To generalize the 
baryons to the large N c case, we follow the prescription suggested by Dashen and Manohar 
H such that there is no valence strange quark in a proton. 

The identity between the large N c QCD and large N c QQCD imposes some constraints on 
the low energy effective field theories of QCD and QQCD, which are the chiral perturbation 
theory (%PT) [pHT2| and the quenched chiral perturbation theory (Q%PT) [Q, respectively. 
%PT and QxPT have the same symmetries and symmetry breaking patterns as QCD and 
QQCD. In both cases, low energy observables are expanded in a power series of the small 
expansion parameters (p, Mqb)/ h. xi where p is the characteristic momentum transfer in the 
problem, Mqb is the goldstone boson mass and A x is an induced chiral perturbation scale. 
In QCD, A x ~ 1 GeV and the goldstone boson masses are M n ~ 140 MeV, M K ~ 500 MeV 




N c — > oo, g — > 0, N c g 2 — > constant, 



(1) 




and M v ~ 550 MeV. In the QCD case (N c = 3 ), xPT gives go 
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where m q is the mass of the flavor q(= u, d, s) quark and we have used = m u . The pion 
decay constant f n = 93 MeV and the renormalization scale \x dependence is absorbed by the 
low energy constants (or counterterms) K 3 and K± defined in | PL Q|| - When N c becomes large, 
f n = 0{\fW c ) ) K4 = 0(1/N C ), while B, quark masses and K 3 are 0(1). Thus as N c — > 00, 
only the B and K 3 terms contribute. 

M*^2m u B{l + 2m u K 3 )+0{m 3 q ). (3) 

Note that the /1 dependence of K 3 is 0(1/N C ), so K 3 becomes /1 independent as N c — > 00. 
In the real world, the chiral logarithms could be numerically enhanced even though they are 
formally higher order in 1/N C . 

In the QQCD (iV c = 3) case, QxPT gives a similar result for the pion mass 

Mf = 2m u B^ 1 1 - ^£joy 2 M 2 log (^f\ + 2m u K? J + • - - , (4) 

where the low energy constants B®, f® and have different values than those in QCD 
and where Mq is the strength of the 7/7/ coupling. In QCD, Mq can be iterated to all orders 
to develop a shifted rj' pole such that rj' becomes non-degenerate with the goldstone bosons 
in the chiral limit. However, in QQCD, all the sea quark loops are dropped such that the 
iteration stops at one insertion of Mq. This gives the v( propagator a double pole dependence 
and makes physical quantities in QQCD more singular than those in QCD in the infrared 
region ||. The quenched pion mass in eq.([|) is an explicit example. The quenched chiral 
logarithm (the Mq term) is more singular than the QCD chiral logarithm in eq.(^). As N c 
becomes large, Mq scales as 1/N C while the other quenched low energy constants scale the 
same way as the corresponding unquenched ones, = 0{VK), and = 0(1). Thus 
as N c — > 00, 

Mf -> 2m u B Q + 2m u K^ + 0(m 3 q ). (5) 

Since QCD and QQCD have the same large -/V c limits, eqs.@ and (§) imply 

B = B Q , K 3 = as iV c -> 00. (6) 

Thus some low energy constants (to be more precise, those that are leading in the 1/N C ex- 
pansion) of xPT and Q^PT are identical in the large N c limit. This implies their differences 
are higher order in 1/N C : 

bQ = b(i + <D(±-) 

K$ =K 3 (l + O(^)) . (7) 

Now we can have a clear idea about how large the quenching error is. The error could 
come from the quenched logarithms which might be formally higher order in the 1/N C 
expansion but numerically significant. Fortunately, those chiral logarithms can be calculated 
using effective field theories so that one can just replace the quenched logarithms by the real 
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ones. Another source of error comes from the difference between the low energy constants 
which is an 0(1/N C ) ~ 30% effect. This is consistent with the ~ 25% quenching error seen 
in simulations (see M, for example). 

Recently kaon weak matrix elements relevant to the K — > 7T7T process have been simulated 
using the quenched approximation by two groups In both calculations, QxPT low 

energy constants are extracted from K —>■ it and K — > vaccum amplitudes in order to 
recover the desired K — > itn amplitude. Based on the above discussion that leads to eq.(^), 
the errors of these analyses are 0(1/N C ) and naively are ~ 30%. Indeed this is the case 
for the AI = 1/2 rule. The simulations give 25.3 ± 1.8 [|5| and 10 ± 4 |3| compared to 
the experimental value 22.2. However, for the case of CP violation e'/e parameter, Ref. || 
obtained (—4.0 ±2.3) x 10~ 4 , consistent with Ref. ||, compared to the current experimental 
average of (17.2 ±1.8) x 10~ 4 . In the e'/e large accidental cancellation happens 

between the 1 = and 1 = 2 contributions such that the estimated ~ 30% errors from 
the two pieces add up to ~ 8 x 10~ 4 . Besides, in the e'/e quenched calculation, there is 



an additional quenched artifact that can be fixed by dropping some "eye diagrams" p3 
Doing this will increase the e'/e by ~ 4 x 10~ 4 |[L4| . Thus in this case the quenching error 
is significant. 

The above examples suggest that in principle, the quenched approximation is typically 
expected to give a ~ 30% error result if there are no big cancellations between different con- 
tributions and if the quenched chiral logarithms (or non-analytic contributions, in general) 
are replaced by the unquenched ones. Now we will try to identify some matrix elements 
that can be reliably computed in the quenched approximation with 100% errors and are still 
worth computing even with that size of uncertainties. 

The first quantity we look at is the proton strange magnetic moment fi s . It is measured in 
the electron-proton and electron-deuteron parity violation experiments The measured 
value has a large uncertainty 

Hf* = [0.01 ± 0.29(stat) ± 0.31 (sys) ± 0.07(theor.)] n.m. (8) 

This implies that the strange quark only contributes around —0.1 ± 5.1% of the proton's 



magnetic moment. In xPT, the first two orders in chiral expansion for fi s is (16| 



^. = A + ^t(^-WF + %F\ (9) 

where D and F are goldstone boson nucleon couplings. Using the standard values for the 
parameters, the second term which is the kaon cloud contribution is ~ 2.0 n.m. This 
combined with the experimental data implies the leading order low energy constant is 
~ —2.0 n.m. Thus a quenched calculation would have a quenching error of ~ 0.8 n.m. 
Again, the large cancellation makes the quenching error relatively big and, hence, the result 
is not reliable within a 100% error. 

The second quantity that is also interesting but even less well measured is the longest 
range parity violating (PV) isovector pion nucleon coupling, h n . For many years, serious at- 
tempts have been made to measure this quantity in many-body systems where the PV effects 
are enhanced. However, currently even its order of magnitude is still not well constrained 



see [[17]] for reviews) largely due to the theoretical uncertainties. Few-nucleon experiments, 
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such as the new np — > d'j experiment at LANSCE fl8|| , are expected to have small the- 
oretical uncertainties and will set a tight constraint on h n . Also, recent studies in single 
nucleon systems in Compton scattering |19| and pion productions |2(J show that the near 
threshold ^fp — > n + n process is both a theoretically clean and experimentally feasible way 
to extract p0| . Given the intensive interests and the level of difficulty involved in the 
experiments, a lattice QCD determination of with a 100% error will be very valuable. In 
fact, issues such as chiral extrapolation has been addressed in []ST| in the context of partially 



quenched approximations |22] in which sea and valence quark masses are different. Since 



one does not expect large cancellations between the first few orders — as what happened in 
the fi s case — we suggest that a quenched simulation of h w is enough to give a 100% error 
result. 

Another interesting PV quantity is the longest range PV photon-nucleon-delta (7iVA) 
coupling. Zhu et al. suggested a model, which was inspired by the hyperon decays, to argue 
that there could be enhanced PV mixing in the initial and final state wave functions such 
that the PV ^yNA coupling is 25-100 times larger than its size from the naive dimensional 
analysis ||23|| . If correct, this enhancement will introduce an easily measurable PV signal in 
^fp —>■ ir + n at the delta peak. To check if the large PV enhancement due to the wave function 
mixing really takes place, one can compute the PV mixing between A(| + ) and iV(1520; § ) 
and between the nucleon(| + ) and iV(1535; \ ) in QQCD. One expects the mixing between 
these states to be 5-20 times larger than their sizes from the naive dimensional analysis 
according to the model of Zhu et al.. This can be checked easily by a quenched calculation 
with a 100% error. 

Recently, Young et al. |24| studied the lattice results of nucleon and delta masses. They 
suggested the difference between the quenched and unquenched cases is largely due to the 
one-loop self-energy contributions calculated in a certain cut-off scheme, while the difference 
due to the low energy constants seems to be smaller than a generic estimation of ~ 30%. 
This might be because in the real simulations, only the mass ratios are computed. The 
absolute mass scale is set by assuming zero quenching error for a certain quantity, then 
the coupling constant and, equivalently, the lattice spacing, is determined accordingly. It 
is conceivable that the ratio of certain quantities can be computed more reliably because 
the leading 0(1 /N c ) errors cancel. An explicit example is the ratio of B meson decay con- 
stants fsj '/b, where B S (B) denotes a low-lying pseudoscalar meson with the same quantum 
number as a bs(bd) state [[|. However, to argue that all the mass ratios have smaller than 
0(1/N C ) quenching errors is beyond the scape of effective field theory. To see if certain error 
suppression really takes place, it would be very interesting to study the N c dependence for 
quenched approximations by direct simulations. 

There is another class of ratios in large N c QCD that the finite N c correction is 0(1/ N^) 
to the leading order. It would be interesting to see whether the correction is still 0(1/ N^) 
in QQCD. 

In summary, quenched QCD and QQCD have the same large N c limit, such that some 
of the low energy constants in their corresponding effective field theories are identical in 
the large N c limit. This implies that some physical quantities can be determined in the 
quenched approximation with 0(1/ N c ) ~ 30% accuracy. A few interesting applications are 
explored. 
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